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We investigate the sea quark polarization in the nucleon by means of a combined next to leading 
order analysis to the recently enlarged set of inclusive and semi-inclusive polarized deep inelastic 
scattering data. Using the Lagrange multiplier method, we asses the uncertainty inherent to the 
extraction of the different spin dependent parton densities in a QCD global fit, and the impact of 
the increased set of semi- inclusive data now available. We comment on future prospects at RHIC 
and JLAB and their potential impact in the future determination of polarized parton densities. 
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For more than fifteen years, polarized inclusive deep inelastic scattering (pDIS) has been the main, if not unique, 
source of information on how the individual partons in the nucleon are polarized at very short distances. Hadronic 
decays, of course, bring us from lower energies hints about the relations between the net polarization of each flavor, 
but only after making strong assumptions about flavor symmetry within the nucleon. These kind of assumptions, 
however, have been seriously challenged by data over the last decade, at least in the realm unpolarized experiments 

una 

The electromagnetic nature of the dominant interaction between the leptonic probe and the target in the pDIS 
experiments performed so far, does not allow to disentangle quark from anti-quark contributions, and thus valence 
quark from sea quarks. Consequently, in spite of strenuous efforts and very successful experimental programs run 
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by several collaborations, our knowledge of the detailed spin dependent parton distribution functions (pPDF) has 
been held hostage to a large extent by our own assumptions. Most notably, the sea quark densities and also that 
I for the gluons. The latter, because it influences the pDIS observables mainly through their energy scale dependence. 
Although at next to leading order accuracy the gluon density contributes to the pDIS cross section directly, this 
O i| contribution is relatively suppressed. 

Many alternative experiments have been conceived, and in some cases already implemented, to improve this situ- 
' ation. The most mature among them are those based on polarized semi- inclusive deep inelastic scattering (pSIDIS), 
, i.e. a pDIS experiment where a particular hadron is tagged in the hadronic final state. Choosing different target 
and final state hadrons, the respective cross sections are sensitive to different combinations of flavored quarks and 
anti-quarks, to be disentangled. These experiments began with the pioneering efforts of SMC, in the late nineties 3], 
followed by those of HERMES at DESY 4], and lately COMPASS at CERN [|, and are planned to be improved at 
the Thomas Jefferson National Laboratory (JLAB) in the near future 

The phenomenological impact of the pSIDIS data proved to be encouraging although scarce in the initial stage: the 
reduced number of data and the relatively large estimated errors, at best allowed to check the consistency between 
pDIS and pSIDIS in a variety of spin- flavor symmetry scenarios 0, @, • With the availability of larger sets of pSIDIS 
data, much more precise, and for final state hadrons and targets of different flavor composition, the situation now has 
changed dramatically. pSIDIS data have a non negligible weight in combined global fits at present, comparable to 
that of inclusive data, and also show clear preferences for the light sea quark polarization. It also helps to constrain 
the strange sea quark and gluon polarization complementing the information already obtained from pDIS. 

Although we can now leave behind the spin-flavor symmetry assumption, translating polarized observables to 
parton distributions always imply the use of some additional piece of information. In the case of pDIS, we need the 
unpolarized structure function F±, or a set of unpolarized parton densities (PDFs), in order to go from the pPDFs 
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to the measured asymmetries. For pSIDIS we also need a set of fragmentation functions (FFs), the functions that 
measure the probability of a quark fragmenting into a given hadron [ToL ITT| , and polarized fracture functions 0] 
if target fragmentation events are taken into account, or a Monte Carlo generator modeling these same processes. 
The results of the global analyses to polarized data show negligible differences when using different sets of modern 
PDFs 0, 0], but the differences are sizable when using different FFs, a feature to be taken into account when 
assessing the uncertainties of the pPDFs. Eventually, pSIDIS data will help to constrain FFs, including it along with 
electron-positron annihilation data in the global fits where they come from. 

In the following we perform a combined next to leading order analysis to the recently updated set of pDIS and 
pSIDIS data, and present the resulting pPDFs. Specifically, we focus on the extraction of sea quark and gluon 
densities, analyzing the constraining power of either set of data on the individual densities. As result, we find not 
only a complete agreement between pDIS and pSIDIS data, but a very useful complementarity, leading to rather well 
constrained densities. 

Using the Lagrange multiplier approach |3fjj|. we explore the profile of the \ 2 function against different degrees 
of polarization in each parton flavor. In this way we obtain estimates for the uncertainty in the net polarization 
of each flavor, and in the parameters of the pPDFs. We compare results obtained with the two most recent sets 
of fragmentation functions, and also within the leading (LO) and next to leading order (NLO) approximation. The 
differences found using alternative sets of FF are found to be within conservative estimates for the uncertainties. 
Nevertheless, there is a clear preference for a given set of FF over the other, shown in a difference of several units in 
the x 2 of the respective global fits. In NLO global fits the overall agreement between theory and the full set of data 
is sensibly higher than in LO case. 

Finally, we analyze the behavior of the cross section for longitudinally polarized proton-proton collisions into neutral 
pions with a wide range of pPDFs sets coming from a rather conservative uncertainty interval. This observable is 
found to be crucially sensitive to the polarized gluon density and therefore an invaluable tool. We compute the 
required precision to be reached in the programed experiments in order to constrain even further this distribution 
and also future sets of pPDFs. A similar analysis is made for forthcoming pSIDIS data to be obtained at JLAB. 



II. CONVENTIONS AND DATA SETS 



Throughout the present analysis, we follow the same conventions and definitions for the polarized inclusive asym- 
metries and parton densities adopted in references 0, 00- however we use more recent inputs, such as unpolarized 
parton densities |l2j | and the respective values for a s . In the totally inclusive case, the spin dependent asymmetries 
are given by Q 



F?{x,Q 2 ) 



<?f(*,Q 2 ) 



F 2 N (x,Q 2 )/2x[l 



R N {x,Q 2 ) 



(1) 



where the inclusive spin-dependent nucleon structure function g^(x,Q 2 ) can be written at NLO as a convolution 
between polarized parton densities for quarks and gluons, Aqi(x,Q 2 ) and Ag(x,Q 2 ), respectively, and coefficient 
functions ACi(cc)P 
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2 ^ q 



A q (x,Q 2 ) + ^fl [ 1 ^\AC q (z)A q (*,Q 2 )+ AC g (z)Ag(*,Q 2 ) 



(2) 



A more detailed discussion about these coefficient functions and their factorization scheme dependence can be found 
in Ref.0|. F^(x,Q 2 ) is the unpolarized nucleon structure function that can be written in terms of F 2 N (x,Q 2 ) and 
R, the ratio of the longitudinal to transverse cross section . The use of the QCD NLO approximations both in 
F^(x, Q 2 ) and in gi(x, Q 2 ) has been shown to reduce the sensitivity of PDFs to higher twists [l7j . 



Analogously, for the semi-inclusive asymmetries we have: 



A? h (x,Q*) U 



J z dzg^ h {x,z,Q 2 ) 
J z dzF^ix^z^Q 2 )' 



(3) 



where the superscript h denotes the hadron detected in the final state, and the variable z is given by the ratio between 
the hadron energy and that of the spectators in the target. The region Z, over which z is integrated, is determined 
by kinematical cuts applied when measuring the asymmetries. 
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For the spin dependent structure function gf f (x,Q 2 ), we use the NLO expression 

g? h (x,z,Q) = \Y& Ui(^Q 2 )d^(^Q 2 ) + ^ 1 ff f§ (A g (|,Q 2 )AC«(M,Q 2 )^?(f,Q 

q,q L x z \ 

+Aq (~ Q 2 ) ACW (x, z, Q 2 )D% (|, Q 2 ) + A 5 (|, Q 2 ) AC« (z, £, Q 2 )Df (J, Q 2 ) | , (4) 

and in order to avoid the convolution integrals in x and z we switch to moment space in both variables as suggested 
in [l9l | and already implemented in [i| . 

For u and d quarks plus anti-quarks densities at the initial scale Qq = 0.5 GeV 2 we propose 

. s Ar x a «(l - xf"(l + j q x s ") 

^ + = N « B(a q + 1 A + 1) + 7, B(a q l 6 q + l,0 q + l) ^ ^ 

where B(a,0) is the standard beta function, while for strange quarks plus anti-quarks we use 

x(As + As) = 2N { — > (6) 
B(a s + l,p s + 1) 

with a similar parametric form for gluons 

The first moments of the quark densities 5q (N q ) 

Sq= [ dxAq (8) 
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are often related to the hyperon beta decay constants F and D through the SU(3) symmetry relations 

6u + Su- Sd- 6d = N u - N d = F + D = 1.2573 (9) 

Su + Su + 5d + Sd - 2(Ss + 5s) = N u + N d - 4N S = 3F - D = 0.579. (10) 

Under such an assumption, the previous equations would strongly constrain the normalization of the quark densities. 
However, as we are not interested in forcing flavor symmetry, we leave aside that strong assumption and relax the 
symmetry relations introducing two parameters, 6b j and esu(3) respectively. These parameters account quantitatively 
for eventual departures from flavor symmetry considerations, including also some uncertainties on the low- a: behavior, 
and higher order corrections, 

N u -N d =(F + D)(l+e Bj ) N u +N d -4N s = (3F-D)(l+e su{3) ), (11) 

and we take them as a measure of the degree of fulfillment of the Bjorken sum rule [2(j and the SU (3) symmetry. 

Equations i|ll[l allow to write the normalization of the three quark flavors in terms of N s , tBj, and £su(3)- Notice 
that no constraints have been imposed on the breaking parameters since we expect them to be fixed by data. The 
remaining parameters are constrained in such a way that positivity with respect to MRST02 parton distributions 
is fulfilled. These last parameterizations are used in order to compute the denominators of equations (1) and (3). 
Similar results are obtained with other modern sets of PDFs. Consistently with the choice for the unpolarized parton 
distributions, in order to compute a s at LO and NLO we use the values of Aqcd obtained Ref.[l^| . 

As anti-quark densities we take 



x a «(l - x)^ 
! %+ !,/%+!) 



xAq = A/g — \ - J ■ rs q = u,d (12) 



for u and d quarks, and we assume s = s since the possibility of discrimination in the s sector is beyond the precision 
of the data (as in the unpolarized case). 
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TABLE I: Inclusive and semi-inclusive data used in the fit. 



Collaboration 


i arget 


Final state 


^ points 


xteis. 


MO 


proton 


inclusive 


10 


Tool 

23] 


SMC 


proton, deuteron 


inclusive 


12, 12 


24] 


E-143 


proton, deuteron 


inclusive 


82, 82 


25] 


[/ i r r 
Ji-lOO 


proton, deuteron 


inclusive 


O A O A 

24, 24 


26] 


Hermes 


proton, deuteron, helium 


inclusive 


o o o 

y, y, y 


Ml 
I4| 


T7 1 A O 


helium 


inclusive 


Q 

8 


27] 


E-154 


helium 


inclusive 


17 


m 


Hall A 


helium 


inclusive 


3 


f29] 


COMPASS 


deuteron 


inclusive 


12 


m 


SMC 


proton,deuteron 


h+, h~ 


24, 24 


® 


Hermes 


proton, deuteron, helium 


h+, h~, 7T+, 7T", K+, K~, K T 


36,63,18 






Total 




478 





Fragmentation functions are taken from either or ^1 > respectively. We also use the flavor symmetry and flavor 
separation criteria proposed in & f the respective initial scales Q 2 . 

DZ + (z, Q 2 ) = Df (z, Q 2 ) = Df (z, Q 2 ) = Df (z, Q 2 ) 

Df (z, Q 2 ) = Df (z, Q 2 ) = Df (z, Ql) = Df (z, Q 2 ) = Df (z, Q 2 ) = (1 - z) Df (z, Q 2 ) (13) 
for the fragmentation into pions, which have shown to be in agreement with unpolarized SIDIS data pH , and 

+ (z, Q 2 ) = D«- (z, Q 2 ) = (l-z) Df (z, Q 2 ) 
Df (z, Ql) = Df (z, Ql) = Df (z, Q?) = Df (z, Q?) = (1 - z) 2 Df (z, Q 2 ) (14) 

for kaons. 

The data sets analyzed include only points with Q 2 > 1 GeV 2 , listed in Table |U and totaling 137, 139, and 37 
points, from proton, deuteron, and helium targets respectively, from pDIS plus 60, 87, and 18, from proton, deuteron, 
and helium targets respectively from pSIDIS. Notice that in the case of HERMES data, we have taken those coming 
from the most recent release, which include corrections for instrumental smearing and radiative effects 



III. GLOBAL ANALYSES 

As it is well known, there are various alternatives for defining the function to be minimized in the global fit [30l l3l| 
To begin with, we consider the most simple and commonly used in fits to polarized data, namely 



where Ei is the measured value of a given observable, Ti is the corresponding theoretical estimate computed with 
a given set of parameters for the pPDFs, and Oi is the error associated with the measurement, usually taken to be 
the addition of the reported statistical and systematic errors in quadrature. This definition ignores the correlations 
between data points from the same measurements, nevertheless it is widely used since in many cases the full correlation 
matrices are not available. In the cases where the correlation matrices between inclusive and semi-inclusive data are 
available, we have taken them into account by means of the appropriate definition for x 2 0. 

In Table ITll we summarize the results of the best NLO and LO global fits to all the data listed in Table[I](478 data 
points). We present fits obtained using alternatively fragmentation functions from reference labeled as KRE, 
and from reference labeled as KKP. Since the fit involves 20 parameters, the number of degrees of freedom for 
these fits is 478-20=458. Consequently, the x 2 values obtained are excellent for NLO fits and very good for LO. The 
better agreement between theory and experiment found at NLO, highlights the importance of the corresponding QCD 
corrections, for the present level of accuracy achieved by the data. 



5 



TABLE II: \ 2 values and first moments for distributions at Q 2 — 10 GeV 2 



set x 2 


2 

Xdis 


2 

XSIDIS 


8u v 


5d v 


5u 


3d 


8s 


Sg 


<5E 


KRE 430.91 
NL ° KKP 436.17 


206.01 


224.90 


0.936 


-0.344 


-0.0487 


-0.0545 


-0.0508 


0.680 


0.284 


205.66 


230.51 


0.700 


-0.255 


0.0866 


-0.107 


-0.0454 


0.574 


0.311 


KRE 457.54 

LO 

KKP 448.71 


213.48 


244.06 


0.697 


-0.248 


-0.0136 


-0.0432 


-0.0415 


0.121 


0.252 


219.72 


228.99 


0.555 


-0.188 


0.0497 


-0.0608 


-0.0365 


0.187 


0.271 



In NLO fits there seems to be better agreement when using KRE fragmentation functions, whereas at LO the 
situation is the opposite. The difference between the total \ 2 values between KRE and KKP NLO fits comes mainly 
from the contributions related to pSIDIS data, while those associated to inclusive data are almost the same, as one 
should expect in a fully consistent scenario. At variance with the NLO situation, in LO fits, the effect of pSIDIS data 
in the fit to inclusive data is strongly dependent on the set FFs used. This suggest a flaw in either the LO description 
of pSIDIS, the LO FFs used, or most likely in both of these ingredients. Notice also that KKP LO fit achieves a \ 2 to 
pSIDIS data close to those obtained at NLO, but the degree of agreement with inclusive data is several units poorer. 
KRE LO fit improves the agreement with inclusive data, but with a higher \ 2 to pSIDIS data. 

Table ITT1 includes also the first moment of each flavor distribution at Q 2 = 10 GeV 2 , and that for the singlet 
distribution <5£, as reference. Most noticeably, while the KRE NLO fit favors the idea of a SU(3) symmetric sea, KKP 
NLO finds u polarized opposite to d and to s. Gluon and strange sea quark polarization are similar in both fits and 
the total polarization carried by quarks is found to be around 30%. At LO, KRE fits also prefers sea quarks polarized 
in the same direction although u is much less polarized than d and s. KKP LO shows u polarized opposite to d and 
to s as in NLO, and the total polarization is around 25%. 

Notice that Au + Au and Ad + Ad, could in principle be determined using inclusive data alone, with no dependence 
on a FFs. However, since pSIDIS data determine both Aq + Aq and Aq, in the combined fit Aq + Aq ends with a small 
dependence on the FFs. The full set of parameters for the different sets can be found in Table Al, in the Appendix. 
A Fortran subroutine providing the scale dependent pPDFs can be obtained upon request from the authors. 




FIG. 1: Inclusive asymmetries computed with LO and NLO pPDFs against the corresponding data. 



In Figures ^ and |2 we show the inclusive and semi- inclusive asymmetries computed with the different parameteriza- 
tions both at LO and NLO accuracy, against the corresponding data sets. The differences between the various sets can 



6 




FIG. 2: Semi-inclusive asymmetries computed with LO and NLO pPDFs against the corresponding data. 

hardly be noticed in the comparison to inclusive data in Figure ^ although are more significant when comparing to 
pSIDIS data, specially in the case of proton targets. This is due to the fact that the main difference between the sets 
are the light sea quark densities, which are probed by pSIDIS processes of proton targets. The pSIDIS asymmetries 
for deuterium targets are, of course, less sensitive to these differences since they average u and d contributions. 

IV. UNCERTAINTIES 

A crucial issue to be addressed before going any further in the interpretation of the results of the previous section 
is to estimate the uncertainties in the extraction of the individual pPDFs by means of the global fit, and also the 
uncertainty that will have any observable computed with them. This has been thoroughly studied in the context of 
unpolarizcd parton distributions, where the number and precision of the data available is much more significant, and 
the correct estimate of the uncertainties arising from PDF in predictions for observables related to new physics is 
mandatory |12ll30l l32| . 

The sources of uncertainty in PDFs are often classified into those associated with experimental errors on the 
data, and those associated with rather theoretical or phenomenological assumptions in the global fitting procedure, 
including: higher order QCD effects in the analyzed cross sections and their scale dependence, the particular choice 
of the parametric form of the distributions at the initial scale, nuclear target corrections, hadronization mechanisms 
and model assumptions such as s = s. Clearly, while the first category is usually under control, the second one is 
particularly difficult to determine. 

Many strategies have been implemented in order to assess the uncertainties in PDFs and their propagation to 
observables, specially those associated with experimental errors in the data. These include the Hessian approach, 
which assumes that the deviation in \ 2 for the global fit is quadratic in the deviations of the parameters specifying the 
input parton distributions, and then propagates linearly these uncertainties from PDFs to observables. Alternatively, 
the Lagrange multiplier method 30] probes the uncertainty in any observable or quantity of interest much more 
directly. This last method relates the range of variation of one or more physical observables dependent upon PDFs to 
the variation in the x 2 used to judge the goodness of the fit to data. Specifically, it can be implemented minimizing 
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the function 

$(\ i ,a j )=x 2 (a j )+Y,*iO i (a j ) (16) 

i 

with respect to the set of parameters aj of the starting PDFs, for fixed values of the Lagrange multipliers A^. Each one 
of the parameters Aj is related to an observable Oi , which is computed from the set of parameters aj . For A; = we 
get the best standard global fit, for which x 2 ( a j) = Xo ancl Oi( a j) = • Varying Aj, and minimizing <I>(Ai, aj), the fit 
to data deteriorates increasing x 2 ( a j) from its minimum while Oi(aj) variates due to the different set of parameters 
aj found. Performing a series of global fits for different values of Xi, we get a profile of X 2 ( a j) f° r a range of values of 
the quantities Oi. In other words, this tell us how much the fit to data deteriorates if we force the PDFs to yield a 
prediction for an observable different to the one obtained with the best fit 0°. 

In order to illustrate the method and apply it to fits of polarized data, in Figure [3] we show the outcome of varying 
the x 2 °f the NLO fits to data against the first moment of the respective polarized parton densities 5q at Q 2 = 10 
GeV 2 , one at a time. This is, to minimize 

= X 2 {aj) +\ q 6q(a,j) q = u,u,d,d, s, g. (17) 

Of course, beyond LO, these first moments are not, strictly speaking, physical observables, however are perfectly well 
defined quantities once the factorization prescription is fixed. The solid lines in Figure |3| correspond to the KRE NLO 
fit (using as input for fragmentation functions those of reference pl|l while the dashed lines correspond KKP NLO 
(using M). 




In and ideal situation, i.e. to have reliable estimates of every source of uncertainty, of correlated experimental and 
theoretical errors, and a quadratic dependence of x 2 m the parameters of the fit, the profile of x 2 would be just a 
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parabola and the la uncertainty in any observable would correspond to A\ 2 = 1 . In order to account for unexpected 
sources of uncertainty, in modern unpolarized global analysis it is customary to consider instead of Ax 2 — 1 between 
a 2% and a 5% variation in x 2 as conservative estimates of the range of uncertainty. 

As expected in the ideal framework, the dependence of x 2 on the first moments of u and d resemble a parabola 
(Figures 3a and 3b). The KKP curves are shifted upward almost six units relative to those from KRE, due to the 
difference in x 2 of their respective best fits. Although this means that the overall goodness of KKP fit is poorer than 
KRE, 8d and 5u seem to be more tightly constrained. The estimates for Sd computed with the respective best fits 
are close and within the Ax 2 = 1 range, suggesting something close to the ideal situation. However for 8u, they only 
overlap allowing a variation in Ax 2 of the order of a 2%. This is a very good example of how the A% 2 = 1 does not 
seem to apply due to an unaccounted source of uncertainty: the differences between the available sets of fragmentation 
functions. 
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FIG. 4: Parton densities at Q 2 = 10 GeV 2 , and the uncertainty bands corresponding to A\ 2 = 1 and A\ 2 = 2% 

An interesting thing to notice is that almost all the variation in x 2 comes from the comparison to pSIDIS data. 
The partial x 2 value computed only with inclusive data, XpDis^ ^ s almost flat reflecting the fact the pDIS data are 
not sensitive to u and d distributions. In Figure 03 we plot XpDis with an offset of 206 units as a dashed-dotted line. 

The situation however changes dramatically when considering Ss or Sg as shown in Figures 3c and 3f, respectively. 
In the case of the variation with respect to 8s, the profile of x 2 is not at all quadratic, and the distribution is much 
more tightly constrained (notice that the scale used for 6s is almost four times smaller than the one used for light 
sea quarks moments) . The XpDis corresponding to inclusive data is more or less indifferent within an interval around 
the best fit value and increases rapidly on the boundaries. This steep increase is related to a positivity constraints on 
As and Ag. pSIDIS data have a similar effect but also helps to define a minimum within the interval. The preferred 
values for 8s obtained from both NLO fits are very close, and in the case of KRE fits, it is also very close to those 
obtained for Su and 8d suggesting SU(3) symmetry. 
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For the total u and d polarization Su + Su and Sd + Sd, the absolute uncertainties are large as shown in Figures 3d 
and 3e respectively, however their values relative to the best fit results are significantly smaller than for sea quarks. 
For 5g, again XpDis J us t defines an interval, as it has been pointed out in previous analyses 0,0, bounded by its own 
positivity constraint in the upper end, and by those of the sea quarks densities, which grow in the lower end. Between 
both ends, there is either some redundancy of parameters or a need of more sensitive data. Again, semi-inclusive data 
help to define a minimum, and it is very close for both fits. Within the A% 2 = 2% uncertainty, the gluon polarization 
estimate is in agreement with a recent analysis of high px hadron production |33j | 

Clearly, in the case of Ss or Sg, the simple minded use of the Hessian approach does not apply. In the case of the 
light sea, the sensitivity to a smaller and internally consistent data set, presumably put us closer to the ideal situation. 

In order to see the effect of the variation in x 2 on the parton distributions themselves, in Figure 01 we show KRE 
best fit densities together with the uncertainty bands corresponding to A% 2 = 1 (darker band) and A\ 2 — 2% (light 
shaded band). As expected, the relative uncertainties in the total quark densities and those strange quarks are rather 
small. For gluon densities the A\ 2 = 1 band is also small, but the most conservative A\ 2 = 2% estimate is much 
more significative. For light sea quarks the A\ 2 — 1 bands are moderate but the A% 2 = 2% are much more larger. 
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FIG. 5: x 2 profiles for LO fits obtained using Lagrange multipliers at Q 2 = 10 GeV 2 . 

The profiles for \ 2 obtained with LO fits show similar features to those obtained at NLO. In Figure [S] we show the 
LO profiles for sea quark densities. Notice that before reaching any conclusion about the uncertainties obtained in 
LO and their comparison with those found in NLO a few remarks are in order. In the first place, the rather large 
difference between the LO and NLO \ 2 values in fits to the same sets of data implies that the uncertainties due to use 
of one or the other approximation, and their related ingredients, are not properly accounted for. Any mean to include 
these uncertainties would certainly reduce the constraining power of the LO fit an thus increase the uncertainties on 
the pPDFs. Unfortunately there is no evident way to do this and also to relate the criteria used to go from profiles 
to uncertainties in NLO and LO. The second point is related to the fact that beyond the LO PDFs are factorization 
scheme dependent and thus the relation between the LO and NLO PDFs and their uncertainties is even less direct. 

The Lagrange multiplier method allows also to analyze the interplay and consistency between inclusive and semi- 
inclusive data Q we mentioned in the previous section. The idea is to apply a Lagrange multiplier to the pSIDIS 
contribution to the function to be minimized, 

$(A,%) = xIdisM + ^xIsidisM (18) 

and perform global fits for different values of A. For A values lower than unity, the weight of pSIDIS data in the fit 
is artificially reduced and the fit becomes increasingly dominated by pDIS data as A — > 0. The partial contribution 
of pDIS data to % 2 , X 2 DIS decreases up to a saturation point given by the best fit to pDIS data only, while that of 
pSIDIS data, XpSiDis? increases. Using A values larger than unity, the fit becomes dominated by pSIDIS data, until 
it saturates at the best fit to pSIDIS data alone. 

In Figure El we show the results for XpDis vs - Xpsmis m n * s w ^h a wide range of values of A in Eq. ED for KRE 
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FIG. 6: xIdis vs. xIsidis for KRE and KKP NLO and LO fits 



NLO (solid line), KKP NLO (dashed line), KRE LO (dotted line) and KKP LO fits (dotted-dashed line), respectively. 
The standard best fit in each curve (A = 1) is denoted by a star. 

In both NLO curves, the XpDis v ^hies obtained in the best standard fits, are very close to the common saturation 
value, corresponding to a fit without pSIDIS data. This suggests a high degree of compatibility between both data 
sets in a NLO framework, since the inclusion of pSIDIS data does not worsens significantly the agreement with pDIS 
data. The increase in the relative weight of pSIDIS data reduces XpSiDis on ^y m a cou ple of units in the case KKP 
fits, what indicates that the fit is strongly constrained and that the best fit result for XpSiDis 1S c ^ ose to its saturation 
value. However, XpSiDis m the best KRE NLO is several units smaller although with a XpDis similar to that of 
KKP. The saturation value obtained in this case is also much smaller. These last features mean that KRE FFs allow 
a much better and flexible fit to pSIDIS data. 

For LO the situation is quite different: the common saturation value for XpDis i s several units larger that the found 
at NLO, what means that even neglecting pSIDIS data the fits obtained at LO are poorer than those of NLO accuracy. 
The best KRE LO fit value for XpDis 1S dose to the saturation value, but at the expense of a rather large XpsiDis 
value. In the case of KKP LO fit, XpsiDis improves significantly but the departure of XpDis f rom its saturation value 
is very large. 



FUTURE PROSPECTS 



One of the measurements most eagerly awaited by the spin physics community is that of single inclusive large_ 
pion production in longitudinally polarized proton-proton collisions, which is right now being run at BNL RHIC 
The spin dependent asymmetry associated to this kind of process, A\ L is defined, as usual, in terms of the ratio 
between the polarized and the unpolarized cross sections, 



Al L = 



dAc 



o^-k" X 



daPP- 



(19) 



which is strongly dependent on the gluon polarization. Indeed, in this observable polarized gluons show up in the 
cross section in the dominant terms. The NLO QCD corrections to this observable have been computed recently 
[35l l3(| finding both large QCD corrections, which make unavoidable the use of the NLO approach, and also a very 
significant dependence on A<?, features that in principle would help to constrain the amount of gluon polarization in 
much more stringent way than in DIS experiments. 

The data obtained up to now by the PHENIX Collaboration suggest a very small asymmetry, consistent with 
pPDFs sets with a moderate gluon polarization. In the following we apply the Lagrange multiplier method in order 



11 



to explore the range of variation of the estimates for this asymmetry associated to the uncertainty in the present 
extraction of pPDFs. 



455 



450 
445 

x 2 

440 



435 



430 



-0.004-0.002 0.002 0.004 0.006 0.008 0.01 



2 A 7T° 

FIG. 7: Profile of the global \ to data against A LL at pr — 6 GeV 

In Figure [7| we show the range of variation of A^[ L at an intermediate value for pt — 6 GeV obtained with different 
sets of pPDFs against the variation of the \ 2 to pDIS and pSIDIS data for these distributions. The profile of x 2 
defines a well defined range of values for A\ L allowed by present pPDFs. The solid line represents the profile of 
X 2 obtained using KRE FFs, both in the global fit to data and in the computation of the asymmetry. The dashed 
line represents the same but for KKP FFs. The minima corresponding to both profiles are very close suggesting a 
cancellation of the associated uncertainty. Notice that both the extraction of the pPDFs, and also the estimate of 
A^ LL with a given set, rely on a set of FFs. The double dependence of the observable on the set of FFs used may had, 
in principle, even potentiated the disagreement. 

Notice that when the fitting procedure explores alternative sets of pPDFs in order to minimize or maximize the 
value of A\ L , not only the gluon distributions varies, but all the distributions. Since the observable is strongly 
dependent on the gluon polarization, there is a very tight correlation between A^ L and Sg, reflected by the fact that 
almost the same sets that maximize/minimize the asymmetry, maximize/minimize Sg, however the role of sea quark 
polarization is non negligible. 

In Figure [7| we have included also the profiles obtained using Sg instead of A n LL in minimization as a dotted line in 
the case of KRE and a dashed dotted line for KKP. Clearly, the sea quark polarization can conspire in order to yield 
a larger/smaller asymmetry than the one obtained with maximum/minimum gluon fit at a given % 2 , effect which in 
much more apparent at larger values of A^ LL . This feature will have to be taken into account in the future for a very 
precise measurement of the gluon polarization. 

Unfortunately, the data collected so far in the first two runs of the PHENIX detector at RHIC have estimated errors 
much larger than the uncertainties in the values A\ L coming from the fits, however this situation is going to change 
dramatically in the near future. In Figure |S] we plot AJ^ L as a function of pt using the best NLO fits coming from 
KRE and KKP FFs (solid lines and dashes, respectively) and also with KRE variants designed to enhance/reduce 
A\ L with Ax 2 = 2% (dotted-dashed lines). We have also included the expected uncertainties for the next two runs 
centered at KRE best fit estimate [37| . Clearly, the future measurement of the asymmetry would certainly be able to 
constrain even further the gluon and therefore reduce the uncertainties in pPDFs. 

Another source of complementary information to further constrain the extraction of pPDFs and also FFs is the 
experimental program of the E04-113 experiment at JLAB, which propose to measure pion and kaon pSIDIS asymme- 
tries for proton and deuteron targets p . In Figure^ we show the profile of x 2 of the global fits using KRE (solid line) 
and KKP (dashes) FFs against the tt~ pSIDIS asymmetry on proton targets at one of the kinematic configurations of 
E04-113 [xb3 — 0.203, Q 2 — 2.3 GeV 2 < z >— 0.5). The expected uncertainty for the measurement of this asymmetry 
shown in the bottom of the figure, is significantly smaller than those of the previous measurements, and also smaller 









1 










V , 




I 


KRE NLO 

KKP NLO 












Pj=6 GeV 



12 



0.03 




0.02 
0.01 

. 7E° 





-0.01 KRENLO 

- - KRE NLO A% 2 =2% 
KKPNLO 

-0.02 . 



2 4 6 8 10 12 14 

P T (GeV/c) 



FIG. 8: A^ L vs. pr computed with different pPDF sets and the estimated errors expected by PHENIX at RHIC. 
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FIG. 9: Profile of the global \ 2 to data against A™ and proposed to be measured by E04-113 at JLAB. 



than the present uncertainty coming from the fit. In Figure ^jp we show the same as in the previous figure but but 
for negative kaons. This asymmetry has not been measured yet and the difference between the predictions coming 
form KRE (solid line) and KKP (dashes) sets is much larger that the expected uncertainties. The measurement of 
this last observable together with the combined effect of data for different targets (proton and deuteron) and final 
state hadrons (positive and negative pions and kaons) will certainly constrain the fit even further, specifically the sea 
quark densities, and at the same time provide a more stringent test on the quark flavor separation for the FFs used 
in the analysis. 
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VI. CONCLUSIONS 



The availability of an important set of new data on polarized processes together with the appropriate theoretical 
tools required to interpret them, obtained and developed in the last few years, have led the extraction of pPDFs in 
the proton to mature and to become an important source of information which will still keep growing in the near 
future. 

In this paper, we have assessed the feasibility of obtaining pPDFs, with special emphasis on the sea quark densities, 
from a combined NLO QCD analysis of pDIS and pSIDIS data. We have estimated the uncertainties associated to 
the extraction of each parton density, finding a well constrained scenario. pSIDIS data is not only consistent with 
pDIS, but improves the constraining power of the fit for all the distributions, being crucial for the light sea quarks. 

Different choices for the FFs used in the analysis of pSIDIS data lead to differences in the light sea quark pPDFs, 
although these differences are within a conservative uncertainty estimate for them. While KRE FFs favor a SU(3) 
flavor symmetric sea, KKP FFs suggest an SU(3) broken one, however the former led to fits of much better quality. 

The first moment of the gluon distribution at a 10 GeV 2 is found to be close to 0.6, constrained to be smaller than 
0.8 and larger than -0.05 within a conservative A^ 2 = 2% range. The upper constraint comming from the requirement 
of positivity of the gluon density and the lower one due to the QCD evolution which forces sea quark densities to 
saturate their own positivity constraints. 

The overall picture found for the quark densities at a 10 GeV 2 is one in which, within uncertainties, up quarks are 
almost 100% polarized parallel to the proton, down quarks anti-parallel in a similar proportion, and sea quarks have 
a small flavor symmetric negative polarization. 

Two programmed experiments, the one based on the PHENIX detector already running at RHIC, and the E04-113 
experiment at JLab will be able to reduce dramatically the uncertainty in both the gluon and the light sea quark 
densities respectively, the latter providing also an even more stringent test on fragmentation functions. 
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Appendix: Parameters of the fits 

We present here the parameters of the best fits found in our analyses, both at LO and NLO. 



Parameter 


KRE NLO 


KKP NLO 


KRE LO 


KKP LO 




-0.0067 


0.0280 


-0.1994 


-0.2413 


esc/(3) 


-0.0122 


-0.0795 


-0.15 


-0.1 


a u 


1.2024 


1.0908 


0.8376 


1.1325 


Pu 


3.4517 


3.2909 


3.5832 


4.3992 


In 


7.4178 


9.2128 


10.270 


13.287 


Su 


1.0722 


0.9956 


1.4483 


2.2157 


ad 


0.6717 


0.4913 


0.6431 


0.8770 


Pd 


4.7090 


4.3816 


4.9575 


5.4356 


Id 


14.999 


19.999 


14.999 


14.999 


5 d 


1.5666 


1.3891 


1.8938 


2.4361 


N s 


-0.0441 


-0.0398 


-0.0367 


-0.0382 


a s 


3.500 


2.5000 


2.1045 


3.4466 


ft = /fe = /% 


11.741 


9.7131 


9.8861 


14.997 


Nu 


-0.0444 


0.0799 


-0.0122 


0.0453 




2.500 


1.1332 


2.5828 


1.1892 


Nd 


-0.0454 


-0.0971 


-0.0408 


- 0.0557 


&d 


0.9970 


0.8979 


0.73443 


1.0639 


N g 


0.1273 


0.0781 


0.1965 


0.082 


Ug 


2.398 


2.2901 


0.0390 


1.320 




2.1398 


1.5698 


2.4051 


4.485 


XTotal 


430.91 


436.17 


457.52 


448.71 


Xdis 


206.01 


205.66 


213.48 


228.99 


XsiDIS 


224.90 


230.51 


219.72 


228.99 



Table Al: Coefficients for NLO and LO best fits 



